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Abstract 

In circuits containing closed loops the operator for the current is determined by charge conser- 
vation up to an arbitrary divergenceless current. In this work we propose a formula to calculate 
the magnetically active circulating current Iting flowing along a quantum ring connected to biased 
leads. By gedanken experiments can be obtained from the response of the 

gran-canonical energy of the ring to an external magnetic flux. The results agree with those of the 
conventional approach in the case of isolated rings. However, for connected rings /ring cannot be 
obtained linear combination of bond currents. 

PACS numbers: 73.63.-b,72.10.-d,73.63.Rt 
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FIG. 1: Sketch of a ring with three sites in contact with two semiinfinite one-dimensional leads. The site of 
the ring are labeled with Roman letters, while the site of the left /right lead are labeled with negative/positive 
integers. 

In this Letter we show that the theory of quantum transport must 
be extended when dealing with circuits containing closed loops in order to calculate the 
magnetic moment which couples to an external magnetic field. In the quantum theory of 

nn 

transport the operator [2|]p[ 

Jmn ^(j^mncln^'n. ^nmCli^iri) (l) 

is interpreted as the electron current operator between sites m and n connected by a bond 

with hopping integral tmn- Such interpretation naturally stems for the continuity equation 

d_ 
dt' 



'"fim. ^ ^ Jmn (2) 



in which the change in density on site m is seen as the sum of the currents fiowing from 
site m to all connected sites n. In a similar way one obtains the formula for the current 
density in a continuum system. We will call Jmn the bond current operator since it de- 
pends on the operators straddling a bond. In the case of a ring, however, and in general 
for circuits containing loops, the continuity alone cannot uniquely fix the current since one 
remains free to add a divergence-less component. Such circuits are not merely accademic. 
The experimental realization of mesoscopic metallic rings Q] has prompted an extraordi- 
nary research activity on the quantum behavior of electrons and fundamental paradimgs 
like Aharonov-Bohm oscillations [9] and persistent currents [l^ are currently under intense 



investigation. Recent progresses in connecting aromatic mo 



ecules to metallic leads have 



brought the ring-like topology into the nano- world as well 

Below we specialize the discussion to tight-binding rings connected to biased leads for 
the sake of definiteness. The continuum counterpart is affected by the same ambiguity 
and deserves a similar discussion. These systems have been mainly considered to study 
the quantum interference pattern of the total current [12, US, Ml and of the ring bond- 

n n t-^' ' J 

currents |15l. Ilq . Nevertheless, scarce attention has so far been given to the calculation 



of the ring magnetic moment. The current pattern along the ring is a superposition of a 
circulating current Jj-ing which is magnetically active and a laminar one. How to calculate 
/ring is the main contribution of this Letter. 

Let us consider the tight-binding model of Fig. [1] described by the Hamiltonian H 

H = tJ2 4cn, (3) 

(m.n) 

where (...) denotes nearest neighbor sites. Due to the ambiguity discussed above it is not 
granted that the physical current flowing through bonds a — b and a — c is the same as the 
expectation value Jab and Jac of the bond current operator in Eq. ([T]), so we must ask how 
the current pattern could be measured. In a macroscopic ring connected to leads, one can 
get the current in each wire by using an amperometer or by exploring the magnetic field 
around each branch of the circuit and performing the line integral. However, for a quantum 
ring this cannot be done and in principle the coherence between alternative paths (that 
electrons explore a — c — b and a — b simultaneously) defies every bond-related definition of 
the circulating current. 

In the case of an isolated ring with Hamiltonian //i.ing(0); the current I^^^^*^ that can 
couple to a magnetic flux and generate the ring magnetic moment isfl?!. 18] 

risolated _ '^(-^ring(0)) 

Here, the magnetic flux is = aab + a^c + ctca where the a's are the phases of the hopping 
integrals t^n = |^mn|e*^?^, in accordance with the Peierls prescription. Experimentally, one 
could get I^^l^^""^ by measuring the torque acting on the ring in a magnetic field. There is 
no ambiguity, since I^f^^^'^'^ = Jab = Jca is the only physical current. 

In this work we define Irmg for the connected ring by a proper magnetic measurement 
to be performed in situ on the ring itself. We shall see that the Hamiltonian contains 
enough information to compute /ring, since the coupling to an external field via the Peierls 
prescription encodes the necessary information. To this end, we must introduce local force 
measurements and illustrate the idea by electric and magnetic thought experiments in par- 
allel, since the two cases illuminate each other. We wish to show that in both cases we need 
a local probe and a local readout of the result. 

Electrostatic experiment. Suppose a macroscopic circuit is prepared in the eigenstate |\l/) 
of the Hamiltonian H. One can gain information about the charge and polarisability at 



some site m of the system by an intensive measurement, e.g. by measuring a force. As a 
probe, one could use a tiny condenser to set a weak electric field of strength D directed, 
say, along the x axis, right at site m. The on-site field shifts the atom by x and changes 
the site energy accordingly, Sm —>■ £m{x) = £m + Dx, while the ground state becomes 
with 1^(0)) = 1^). Now, we must decide precisely which force to measure as a response to 
the local probe. If the whole circuit could be treated like a rigid body, one could measure 
F — — ^Ij-^o where E{x) — {^{x)\H{x)\^{x)) is the total energy; F is total force on the 
system. By exploiting the Hellmann-Feynman theorem one finds that the observable is 
—F/D = Um = (^|?^m|^)j i-e., the average density on site m, which is often obtainable 
by simpler means. However, as the probing electric field is local, a local readout of the 
experiment is desirable. The on-site force on the atom can be measured using, e.g., an 
atomic force microscope and can be expressed in terms of the local energy as 

^=-^Uo. (5) 

where 

E.t{x) = {^{x)\{em{x) - (6) 

is the grand-canonical energy of site m with chemical potential /x. The external circuit works 

as a reservoir for particles and heat. The use of the grand-canonical formalism ensures the 
gauge invariance of the theory versus shifts of the energy origin. The local force measurement 
yields F — —n^^D where n^* — {^\hm\^) -\- Srim with 

Snm = '-^^ (*(^)l^-l^(^))L=o- (7) 

The extra contribution 6nm is interesting since 6nm = {^m — f^)Xm where Xm = 
h'Sx (^(•^)l^nil^(^))L=o ^^'^ ratio of the polarization charge d{'^ {x)\hm\'^ (x)) to the ex- 
ternal potential Ddx induced by a small shift of the atom in the field. Thus Xm brings 
information on the local dielectric response, while the factor Sj^ — ^ accounts for the work 
done to bring charge from infinity to site m. For example, if the overall charge on the atom 
is negative, the field will shift it towards positive potentials, hence the site will be more 
attractive for electrons and Dx < 0; thus we may predict that d{^ {x)\n„i\^ {x)) > and 
Xm < 0. Then if Em < A*, that is, the level is more than half filled, we may expect a positive 
Srim, and a further increase of the electron population, while less than half filled levels will 



4 



tend to be emptied. In Fig. [2] we display the ground state response Xm as well as the den- 
sities Tim and n^^^ for a one- dimensional tight-binding chain with nearest neighbor hopping 
t and zero on-site energy, Em = 0, as a function of the chemical potential fi. Note that at 
high filling the local response of the system does not vanish due to a split-off state at energy 
larger than 2. 




FIG. 2: (Color online) Local response Xm for a one-dimensional tight-binding chain as a function 
of the chemical potential for Em = 0. The inset shows the ground state local density rim (red, 
dashed) and the total response density n|:^* (blue, solid) as in Eq. ([7]). 

Magnetic experiment. We develop the magnetic case by a gedanken experiment as far as 
possible in parallel with the above, aiming at a definition of the ring current that corresponds 
to Eq. ([7]). While the leads are biased and the current flows, we switch a weak uniform 
magnetic fleld of strength B forming an angle 6 with respect to the normal of the ring. 
We denote with the magnetic flux through the ring and with |^E'(0)) the current carrying 
state of the system after all transient effects have disappeared. The derivative of the total 
energy E{4>) = (^(0)|^(0)|\Ef(0)) would give the torque acting on the whole system. As 
in the electrostatic experiment, the wave function is modifled everywhere even though the 
magnetic perturbation is localized and hence (^^(0)1^(0) — ifi.ing(0)|\E'(0)) depends on the 
flux and the external circuit experiences a torque as well. A further problem in using the 
variation of the total energy is related to the choices of the Peierls phases along the ring. The 
definition of the local torque should only depend on the magnetic flux = aab + cnhc + «ca 
through the ring. However, one can easily reahze that the choice aah = and abc = O-ca = 
(cl) and the choice aab = «bc = and aca = (c2), see Fig. [3l are not related via a gauge 
transformation and hence lead to different derivatives of the total energy with respect to 0. 
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The asymptotic current- carrying eigenstate 16| of H can be used in both cases to invoke the 



Hellmann-Feynman theorem. The choice cl then leads to the bond current Jab while the 
choice c2 to Jca- The use of the total energy to calculate the torque experienced by the ring 
is therefore ambiguous. Next, we show that a local measurement on the ring is much more 
rewarding. 

The ring interacts with the magnetic field and its energy relative to the chemical potential 
11 is 

^ring(0) = (^(0)|^ring(0) -/UiVring|^(0)) 

= -MB COS 9 = -M^, (8) 

where A'^^ng is the number operator of the ring, S is the ring surface and M is the ring 
magnetic moment. We have no information about S, however, we may say that 

1 __M_dE^ 

The definition of the ring current does not suffer from the ambiguity originating from dif- 
ferent possible choices of the Peierls phases. This follows from the fact that the operator 
Hnng{(p) — /Wiring is a local Operator and hence its average only depends on the projection 
onto the ring of the single particle states {ipk} forming the Slater determinant Let 
us discuss this crucial point in more detail. We consider the choice cl and let ipk{f^) be the 
amplitude on site m of the /c-th one-particle eigenstate of the total Hamiltonian, see Fig. 

Similarly we denote with ipki^n) the one-particle eigenstate corresponding to the choice 
c2. The choice c2 can be transformed in a gauge equivalent phase configuration (choice c2') 
in which the Peierls phases along the ring are the same as those in cl but in the right lead 
there is a bond, e.g., the bond 1 — 2, that acquires a phase —0. The eigenstate ipkijn) of 
choice c2 is transformed in the eigenstate ip'i^{m) accordingly. It is straightforward to realize 
that ip't^im) = e^^ipki'tTT') for m > 1 while ip'f.i'm) = ipk{fn) otherwise. As a consequence the 
average of the local operator Hring{4') — fJ'Nrmg over the gauge inequivalent configurations cl 
and c2 does not change. Such result is independent of the choice of the Peiers phases along 
the ring provided that aab + dbc + dca = 0- 

Results and discussion. To calculate the ring current from Eq. we use an embedding 
technique. We consider the system of Fig. [T] with zero on site energies everywhere and 
hopping \t\ between the sites connected by a link. Let Ul/ji be the bias applied to the 
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FIG. 3: Configuration of the Peierls phases along the bonds of the system. For choice cl the phase 
is localized on the bond a — b while for choice c2 the phase is along the bond a — c. Starting from c2, 
the multiplication of the fermion operators on sites c, b, and —1 by e'^'^ is a gauge transformation 
which corresponds to the configuration c2'. 

left/right lead and h be the matrix of the one-body operator Hring{<P) = ^iji'P)clcj. 
Then the average i?i.mg(0) in Eq- dH]) can be expressed in terms of the lesser Green's function 
G< as 

^Tr,i,g[(/i-/x)G<H], (10) 

where the trace is taken over the sites of the ring. The matrix C^lu) is the Fourier transform 
of the lesser Green's function G^{t,t') for times — > oo and can be written as C^luj) = 
G^{uj)[I]^{uj) + The retarded/advanced Green's function projected onto the 

ring is G^/^{u) = [uj - h - S^/A(cj)]-i with S^/^ = T.^^^ + T^^^ the retarded/advanced 
embedding self-energy of the left and right leads. Using the fluctuation-dissipation theorem 
for lead a = L,R one obtains for the lesser embedding self-energy Sa(co') = —2if{u — 
f/a)Im[E^(a;)], where f{uj) is the Fermi distribution function at chemical potential /i. The 
retarded and advanced components are related as = [S^]^ with a = L,R. For one- 
dimensional tight-binding leads with nearest neighbor hopping t the self-energies have only 
one nonvanishing matrix element, namely [S^(co')]jj = Sia^jaCriuj — Ul) and [S^(co')]ij = 
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6ibSjb(j{uj — Ur). The function a{uj) can be easily calculated and reads 



+ IT]) + 2t 



1 + 



At 



(11) 



{u+iri)-2t 

The ring current in Eq.Q is obtained by taking the flux derivative of -E'ring(0) in Eq- (|T0|) in 
= 0. The flux derivative of the one-body matrix hamiltonian h yields a linear combination 
of the one-body matrix bond currents with coefficients oibc/ 'P and aca/<P- This term 

alone would then be dependent on the Peierls phase configuration. The independence of 

reads 



ring from the phase configuration is restored by adding the flux derivative of G'^{uj) which 



d(j) \o(p 



dcp J 



(12) 



Thus /ring can be expressed in terms of ^10=0 and Green's functions at = 0. We wish to 
emphasize that for = the ring current has no diamagnetic contribution. 
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FIG. 4: (Color online) Plot of the bond currents Job (green, dot-dashed) and Jca (red, dashed) and 
of the ring current /ring (black, solid) defined in Eq. ([9]) as a function of the applied bias voltage 
\Jr = —Ul = U for different values of the equilibrium chemical potential /i. Energies are in units 
of the hopping parameter t and currents are in units of IGq with Go the quantum of conductance. 

In Fig. m we display the I/V characteristic of the ring current as well as of the bond 
currents along the a — b and c — a bonds for different values of the chemical potential 
/i = —1.5, 0. — 1.0, 0.5, 1.0 in units of t. The bond-currents are computed using a Landauer- 



like formula derived in Ref. 



id In all cases the ring current is quadratic in U for small U 



meaning that the ring conductance is always zero. 



On the contrary, the bond currents are generally linear in U except for special values of 
fi at which the bond conductance vanishes [iQ]. One of these special values of /i is yU = 1. In 
this case the ring current coincides with the bond current along the c — a bond for all values 
of U. We also observe that /ring has a maximum as a function of U and that the position of 
the maximum shifts towards high bias by increasing the chemical potential. For negative n 
the maximum is located in correspondence of the maximum of the bond currents while for 
positive fi is locate in correspondence of their minimum. All currents correctly vanish for 
bias [/ = 2 as the left continuum is lifted by 2 while the right continuum is lowered by the 
same amount. Since the bandwidth is 4 the bias U = 2 represents the minimum value of U 
for which there is no overlap between the left and right continua. 

In a similar way one can compute /ring for rings with sites, arms of different length 
and different hopping as well as onsite energy parameters. We have verified that /ring = 
for rings symmetrically connected, as physically expected. 

In conclusion we have pointed out that for quantum circuits containing closed loops the 
theory of quantum transport needs to be extended in order to compute the loop magnetic 
moments which couple to an external magnetic field. By a suitable gedanken experiment we 
have been able to define a ring current which is independent of the Peierls phase configuration 
provided that + + «ca = is kept constant. The explicit calculation of /ring in a 
ring connected to one-dimensional tight-binding leads show that /^ng is, in general, not 
given by a linear combination of the bond currents, even though there are common features. 
Our procedure differs from the one of Ref. [l3] where the magnetic moment is obtained 
by an average over the bond currents. We believe that our treatment can be compared 
with experiment by measuring the ring magnetic moments, and paves the way to include 
induction and self-induction effects in quantum transport theory. In the extended theory, 
even in the absence of an external magnetic field one will need to consider a flux = cLI^ng 
where L is the self-induction coefficient. 
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